Collective Excitations in the Charge-Ordered Phase of a-(BEDT-TTF)2l3 
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The charge response of charge-ordered state in the organic conductor a-(BEDT-TTF)2l3 is char- 
acterized by dc resistivity, dielectric and optical spectroscopy in different crystallographic directions 
within the two-dimensional conduction layer. Two dielectric modes are detected. The large mode 
is related to the phason-like excitation of the 2A;f bond-charge density wave which forms in the ab 
plane. The small dielectric mode is associated with the motion of domain-wall pairs along the a- 
and b-axes between two types of domains which are created due to inversion symmetry breaking. 

PACS numbers: 71.27.-|-a, 71.45.Lr, 77.22. Gm, 74.25.N- 
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Electron-electron and clcctron-phonon interactions, in 
particular in the systems with reduced dimensionality, 
are known to be the driving force for the formation of new 
ordered states of matter P, Q • Among the most intrigu- 
ing phenomena found in these systems are broken sym- 
metry phases like charge- and spin-density waves (CDW, 
SDW), charge order (CO), antiferromagnetic and spin- 
Peierls phases; all of them show a large variety of non- 
linear properties and complex dynamics, including collec- 
tive excitations ■ This wealth of charge-ordered phe- 
nomena in one- (ID) and two-dimensional (2D) strongly 
correlated systems is brought in by the variety in lattice 
structure, represented by anisotropic networks, both in 
the electron hopping t and in the inter-site Coulomb in- 
teraction V. While CO is an effect observed both in inor- 
ganic and organic materials with strong electronic corre- 
lations, the organics present a more robust and "clean" 
state, where in some cases CO does not compete with 
other ground states, as opposed to high-temperature su- 
perconducting materials, for instance [7|. 

Theoretical considerations indicate that the charge dis- 
proportionation is driven by Coulomb repulsion. In par- 
ticular, for systems with a quarter-filled conduction band 
even large values of on-site repulsion U are not suffi- 
cient to transform the ground state from metallic to 
insulating; rather the inter-site Coulomb interaction V 
is responsible to stabilize a Wigner-crystal type phase 
(Seo et al. in Ref. Q). Indeed, there are several well- 
established charge-ordered organic ID and 2D systems 
Q : charge modulation driven by long-range Coulomb re- 
pulsion was first evidenced in the ID organic conductors 
pi-DCNQI)2AgbyKanodainRef. % and (TMTTF)2A 
and also in 2D conductors based on the BEDT- 
TTF [bis(ethylenedithio)tetrathiafulvalene] molecule: 9- 
(BEDT-TTF)2RbZn(SCN)4 and a-(BEDT-TTF)2l3 by 
Takahashi et al. [8]. 

The organic conductor a-(BEDT-TTF)2l3 (a-ETzIa) 
is the most prominent example of charge order in 2D or- 
ganic conductors. The crystals are formed by alternating 



anion and donor layers in the ah plane; the unit cell is 
triclinic and contains four BEDT-TTF molecules. The 
BEDT-TTF layer consists of two types of stacks: Stack 
I is weakly dimerized and composed of crystallographi- 
cally equivalent molecules A and A', while the stack II 
is a uniform chain composed of B and C molecules [loj . 
At high temperatures the system is a semimetal with 
small electron and hole pockets at the Fermi surface 11 1 . 
At the metal-to-insulator transition Tco — 136 K the 
conductivity drops by several orders of magnitude and a 
temperature-dependent gap opens in the charge and spin 
sector which indicates the insulating and diamagnetic na- 
ture of t he g round state [13] ■ Nuclear magnetic resonance 
(NMR) |l2!| and synchrotron X-ray diffraction measure- 
ments [l3i] demonstrate that the charge order, whose fluc- 
tuations are already observed at high temperatures [14| . 
develops at extended length scales below Tco- The es- 
timated charge values of the molecules are A = 0.82(9), 
A' = 0.29(9), B = 0.73(9) and C = 0.26(9). The ex- 
act site assignment is not settled yet since these values 
slightly differ from those found in the NMR, vibrational 
infrared and Raman spectroscopy 13, [3 ■ Observed 
structural changes comprise two effects. First is a sym- 
metry reduction at the phase transition with the space 
group PI changing into PI in the CO state, which im- 
plies four non-equivalent BEDT-TTF molecules in the 
unit cell. Second is the modulation of overlap integrals 
due to changes in dihedral angles at low temperatures. 
In other words, the charge order comprises "horizontal" 
charge stripes along the 6 crystallographic axis of charge- 
poor (CP) sites, the A' and C molecules, and charge-rich 
(CR) sites, the A and B molecules, together with a bond 
modulation. The latter infers that the nature of CO in 
this system should not be regarded as fully localized; in- 
stead, a CDW picture would be more suitable. 

The considerable knowledge of the CO pattern in a- 
ET2I3 makes it a suitable candidate in the search for 
collective excitations in this CO phase. Early results 
indicated the existence of a broad relaxation in radio- 
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Figure 1: (a) Frequency dependence of the real (e') (inset) and imaginary (e") (main panel) part of the dielectric function in 
Q:-ET2l3 at representative temperatures for E || [110]. Below 75 K two relaxation modes are observed - full lines for 47 K show 
a fit to a sum of two generalized Debye functions; dashed lines represent contributions of the two modes. Above 75 K only 
one mode is detected, and the full lines represent fits to a single generalized Debye function, (b) Dielectric strength (upper 
panel) and mean relaxation time with dc resistivity (points and line, respectively, lower panel) in a-ET2l3 as a function of 
inverse temperature, for E || [110]. (c) Anisotropy of large dielectric mode mean relaxation times in a-ET2l3 (points, left 
axis). The temperature behavior closely follows the dc conductivity anisotropy (solid line, right axis). In the inset the mean 
relaxation time is plotted as a function of inverse temperature; full and empty symbols mark the large and small dielectric 
mode, respectively, for E along the a- (circles) and 6-axis (triangles). 



frequency range with a large dielectric constant of the 
order of 10^ as well as sample-dependent nonlinearities 
[l7| . In an attempt to characterize the collective ex- 
citations of the charge order, we have undertaken dc 
and ac conductivity- anisotropy measurements on care- 
fully oriented single crystals of Q;-ET2l3. We discovered a 
complex and anisotropic dispersion in the charge-ordered 
state of a 2D organic crystal. First, similar to the Peierls 
CDW state, we observe long-wavelength charge excita- 
tions with an anisotropic phason-like dispersion, which 
we detect as broad screened relaxation modes along both 
the a- and 6-axis of the 2D BEDT-TTF plane. Sec- 
ond, we detect short wavelength charge excitations in the 
form of domain wall pairs, created by breaking the inver- 
sion symmetry, which are less mobile and induce a much 
weaker polarization, again along both crystallographic 
axes. Also, we observe both types of excitations along 
diagonal direction of BEDT-TTF plane. Our results are 
well understood within the theoretical model by Clay et 



al. [18[ which shows that the CO phase with horizontal 



stripes of localized charges is best described in terms of 
a bond-charge density wave. 

The dc resistivity of Q;-ET2l3 was measured between 
room temperature and 40 K. In the frequency range 
0.01 Hz-lO MHz the complex dielectric function was ob- 
tained from the complex conductance measured by two 
setups. At high frequencies (40 Hz-10 MHz) an Agilent 
4294A precision impedance analyser was used. At low 
frequencies (0.01 Hz-3 kHz) a setup for measuring high- 
impedance samples was used based on lock-in technique. 
At frequencies 6-10000 cm~^ the complex dielectric func- 
tion was obtained by a Kramers-Kronig analysis of the 



infrared reflectivity. All experiments were done on flat, 
planar high-quality single crystals along three directions 
within the plane: in the a and b crystallographic axes as 
well as in [lIO] (diagonal) direction To!]. An influence 
of extrinsic effects in the dielectric measurements, espe- 
cially those due to contact resistance and surface layer 
capacitance, was ruled out with scrutiny [20| . 

Fig.Ufa) shows representative spectra at three selected 
temperatures for E || [110] taken below Tqo — 136.2 K 
in the charge-ordered state. Most notably, between 35 K 
and up to 75 K two dielectric relaxation modes are dis- 
cerned. The complex dielectric spectra e(w) can be de- 
scribed by the sum of two generalized Debye functions 
eiLj) - eoo = AeLD/[l + (iwro.Lo)'""^"] + A£sd/[1 + 
(jajTo,sD)^~"^°] where £oo is the high-frequency dielectric 
constant, Ae is the dielectric strength, tq the mean re- 
laxation time and I — a the symmetric broadening of the 
relaxation time distribution function of the large (LD) 
and small (SD) dielectric mode. The broadening param- 
eter 1 — a of both modes is typically 0.70 ± 0.05. In Fig. 
[Ijb) the dielectric strengths and mean relaxation times 
are plotted as a function of inverse temperature. The 
dielectric strength of both modes does not change signif- 
icantly with temperature (Aeld ~ 5000, Aeso ~ 400). 
At approximately 75 K the large dielectric mode overlaps 
the small mode. It is not clear whether the small dielec- 
tric mode disappears at this temperature or is merely 
obscured by the large dielectric mode due to its relative 
size. However, above 100 K, when the large dielectric 
mode shifts to high enough frequencies, no indication can 
be found of a smaller mode centered in the range 10^- 
10^ Hz. Accordingly, above 75 K fits to only one gen- 
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eralized Debye function are performed which we identiiy 
with the continuation of the large dielectric mode. All 
three of the large-mode parameters can be extracted in 
full detail until it exits our frequency window at about 
130 K. At temperatures up to 135 K (just below Too) 
we can determine only the dielectric relaxation strength 
by measuring the capacitance at 1 MHz. 

The most intriguing result is that the dielectric 
strength and temperature behavior of the mean relax- 
ation times differ greatly between the two dielectric 
modes. The large one follows a thermally activated be- 
havior similar to the dc resistivity, whereas the small 
mode is almost temperature-independent. For both ori- 
entation E II a and E || 6 the results are comparable 
to the findings along the E || [110]. There is no pro- 
nounced anisotropy or temperature dependence in the 
dielectric strength: the Ae values of both the large and 
small dielectric modes (not shown) correspond to those 
measured with E || [110]. However, an anisotropy in 
tq.ld is clearly visible and its evolution closely follows 
the dc conductivity anisotropy [Fig.[TIc)]. A similar con- 
ductivity anisotropy has been observed in the CO phase 
of (TMTTF)2AsF6 [llj. In our samples, despite the 



temperature-dependent activation, the anisotropic trans- 
port gap in the CO phase for E || a and E || 6 can be 
estimated to about 2A = 80 meV and 40 meV. Con- 
versely, our optical measurements for E a and E || 
reveal the isotropic gap of about 75 meV [22]. This is not 
surprising since systems with complex band structure can 
exhibit distinct optical and transport gaps: optical mea- 
surements probe direct transitions between the valence 
and conduction band, while dc transport is governed by 
transitions with the smallest energy difference between 
the two bands. 

The observed ac conductivity data demonstrate a com- 
plex and anisotropic dispersion in the CO state. First, 
similar to CDW state, we observe broad screened relax- 
ation (large dielectric) mode along the a- and 6-axis and 
for E II [lIO]. This fact indicates that the observed 
dielectric relaxation is not associated with stripe ori- 
entation and its motion along one of the two in-plane 
crystallographic axes. Rather, this mode can be inter- 
preted as a signature of long-wavelength CDW excita- 
tions possessing an anisotropic phason-like dispersion. A 
2D dispersion with similar dielectric strength amplitude 
and relaxation time behavior was previously found in 
the commensurate CDW phase developed in ladder lay- 



ers of Sri4Cu2404i [23|-l25|. In Q;-ET2l3 Kakiuchi et al. 



suggested a 2/cf CDW formation along the zig-zag path 
CABA'C of large overlap integrals detected in their X- 



ray diffraction measurements [13|. A theoretical model 
for the related quarter-filled 0-ET2X systems may pro- 
vide additional insight Clay et al. showed that the 
CO phase with horizontal stripe phase is characterized by 
a 1100 modulation of site charges along the two indepen- 




Figure 2: Schematic representation of a 2A:f bond-CDW in 
a-ET2l3. Thin dotted line denotes a unit cell. Single and 
double lines represent modulated overlap integrals along two 
p-directions (see text). 



and with a 1010 modulation along the b axis perpendicu- 
lar to the stripes. In addition, this CO is accompanied by 
a modulation of the overlap integrals along the same p- 
direction. In other words, such a CO phase corresponds 
to the 2/cf modulation of bonds and site charges, e.g., a 
combined bond-CDW along the two BEDT-TTF plane 
p-directions. An analogous, albeit more complex, 2fcF 
modulation of overlap integrals indeed develops along the 
p-directions of a-ET2l3, ACA'BA and ABA'CA (see Fig. 
H]) [1^ . It is plausible to look for the origin of phason-like 
dielectric relaxation in such a 2A;f bond-CDW. In this 
case the energy scale of barrier heights is close to the 
single-particle activation energy indicating that screen- 
ing by single carriers responsible for the dc transport 
is effective for this relaxation. The fact that the tem- 
perature behavior of tq^ld anisotropy closely follows the 
dc conductivity anisotropy has important implications: 
while the CDW motion is responsible for the dielectric 
response, the single electron/hole motion along the two 
p-directions, possibly zig-zagging between them, is re- 
sponsible for the observed dc charge transport. 

Second, we address the small dielectric mode. Here we 
note the twinned nature of the CO phase due to breaking 
the inversion symmetry and formation of a 1010 modu- 
lation of site charges along the a-axis, with one domain 
being (A,B)-rich and the other (A',B)-rich 13|. Indeed, 
the ferroelectric aspect of the CO phase was also the- 
oretically suggested 



and experimentally probed [26 1. 
Our data can be naturally attributed to the motion of 
charged kink-type defects - solitons or domain walls in 
the CO texture. Charge neutrality constraint of the CO 
in Q;-ET2l3 (a change of stripes equivalent to strictly re- 
placing unit cells of one twin type with another) suggests 
two types of solitons and/or domain walls. The first one 
is the domain wall in pairs (a soliton-antisoliton pair) be- 
tween CR and CP stripes along the 6-axis, which we get 
if we impose the constraint along the 5-axis [Fig. [SJa)]. 
The second type of domain-wall pair is given by applying 
the constraint along the a-axis in such a way that the do- 
main walls' interior contains both charge signs [Fig.|31[b)]. 



dent p-directions parallel to the larger overlap integral. The motion of such entities induces a displacement cur- 
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Figure 3: Two different types of domain wall pairs in the 
charge-ordered phase of a-ET2l3. (A,B)- and (A',B)-rich 
unit cells are symbolically represented as +— or — h cells 
which form CO stripes. For simplicity we omit the B and 
C molecules. Gray thick lines stand for charge-rich stripes. 
Thin black lines denote a domain wall pair. 



rent and can therefore be considered as the microscopic 
origin of polarization in the CO state. Their relaxation, 
being nearly temperature-independent, cannot be dom- 
inated by resistive dissipation, rather it is governed by 
low-energy barriers, similarly as observed in ferromag- 
netic domain phase Finally, the observed dielectric 
strength As of the small mode of about 400 confirms the 
domain- wall assignment since their dynamics, in contrast 
to phason-like one, is commonly found to be character- 
ized by much smaller dielectric constants (of the order 
1000 and less) 

In conclusion, we demonstrate the development of 
in-plane anisotropic conductivity accompanied by the 
appearance of two low-frequency dielectric relaxation 
modes in the charge-ordered phase with a horizontal 
stripe structure of a-(BEDT-TTF)2l3. The large dielec- 
tric mode features an anisotropic phason-like behavior 
which we associate with the 2kp bond-charge density 
wave formed within the two p-directions in the BEDT- 
TTF plane. On the other hand, we ascribe the small 
dielectric mode to the motion of domain-wall pairs at in- 
terfaces between two different types of charge order do- 
mains created due to inversion symmetry breaking. This 
unusual appearance of a phason-like relaxation alongside 
a soliton-like mode underlines the complexity of collec- 
tive excitations in such a charge order system. Our re- 
sults strongly suggest that a localized Wigner-like picture 
in terms of stripes of localized site charges is not appro- 
priate for charge order in Q!-(BEDT-TTF)2l3. Instead, 
our results are consistent with a description in terms of a 
2D bond-charge density wave. Our findings demonstrate 
that the issue of collective excitations in broken symme- 
try phases becomes ever more subtle and call for further 
work in order to clarify whether such a dispersion is com- 



mon to 2D charge order phases with horizontal stripes in 
diverse strongly correlated systems. 

We acknowledge S. Mazumdar for valuable and illu- 
minating discussions. We thank G. Untereiner for the 
sample preparation and T. Vuletic for his help in data 
analysis. This work was supported by the Croatian Min- 
istry of Science, Education and Sports under grants 035- 
0000000-2836 and 035-0352843-2844 and by the Deutsche 
Forschungsgemeinschaft under grant DR 228/29-1. 



[1 

[2: 

[3: 

[5: 
[6: 
[7: 
[s: 

[9: 
[10: 

[11 

[12 
[13 
[14 
[15 

[16 
[17 

[18 
[19 

[20 
[21 
[22; 

[23' 
[24 
[25^ 
[26' 
[27 
[28' 
[29 



Electronic address: ltivek@ifs.hrl 

URL: http : //real- science . if s .hr/| 

P. Fulde, Electron Correlations in Molecules and Solids, 
S"^"* edition. Springer- Verlag, Berlin (1995). 
Molecular Conductors, Thematic Issue, Chemical Re- 
views 104, No. 11, November 2004. 

G. Griiner, Rev. Mod. Phys. 60, 1129 (1988); 66, 1 
(1994). 

P. B. Littlewood, Phys. Rev. B 36, 3108 (1987). 

R. J. Cava et al. Phys. Rev. B 31, 8325 (1985). 

T. Vuletic et al. Phys. Rep. 428, 169 (2006). 

J. M. Tranquada et al. Phys. Rev. B 78, 174529 (2008). 

Organic Conductors, Special Topics Section, J. Phys. 

Soc. Jpn. 75, No. 5, May 2006. 

P. Monceau, F. Y. Nad, S. Brazovskii, Phys. Rev. Lett. 
86, 4080 (2001). 

K. Bender et al. Mol. Cryst. Liq. Cryst. 107, 45 (1984); 

K. Bender et al. ibid 108, 359 (1984). 

T. Mori et al. Chem. Lett. 1984, 957 (1984). 

Y. Takano et al. J. Phys. Chem. Sohds 62, 393 (2001). 

T. Kakiuchi et al. J. Phys. Soc. Jpn. 76, 113702 (2007). 

S. Moroto et al. J. Phys. IV (France) 114, 339 (2004). 

J. Moldenhauer et al. Synth. Met. 60, 31 (1993); M. Dres- 

sel and N. Drichko, Chem. Rev. 104, 5689 (2004). 

R. Wojciechowski et al. Phys. Rev. B 67, 224105 (2003). 

M. Dressel et al. J. de Physique I (France) 4, 579 (1994); 

M. Dressel et al. Synth. Met. 70, 929 (1995). 

R. T. Clay et al. J. Phys. Soc. Jpn. 71, 1816 (2002). 

Single crystals were oriented in the mid-infrared set-up 

before, and in the X-ray set-up after experiments. 

T. Ivek et al. Phys. Rev. B 78, 035110 (2008). 

B. Korin-Hamzic et al. Phys. Rev. B 73, 115102 (2006). 

C. Clauss, N. Drichko, D. Schweitzer, and M. Dressel, 
Physica B (2009), doi:10.1016/j.physb.2009.11.036. 

T. Vuletic et al. Phys. Rev. Lett 90, 257002 (2003). 

T. Vuletic et al. Phys. Rev. B 71, 012508 (2005). 

P. Abbamonte et al. Nature 431 1078 (2004). 

K. Yamamoto et al. J. Phys. Soc. Jpn. 77, 074709 (2008). 

M. Pinteric et al. Eur. Phys. J. B 11, 217-225 (1999). 

Y. Tokura et al. Phys. Rev. Lett. 63, 2405 (1989). 

S. Tomic et al. Synthetic Metals 120, 695-698 (2001). 



